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The impact of control sequences on the environmental coupling of a quantum system can be 
described in terms of a filter. Here we analyze how the coherent evolution of two interacting spins 
subject to periodic control pulses, at the example of a nitrogen vacancy center coupled to a nuclear 
spin, can be described in the filter framework in both the weak and the strong coupling limit. A 
universal functional dependence around the filter resonances then allows for tuning the coupling 
type and strength. Originally limited to small rotation angles, we show how the validity range of 
the filter description can be extended to the long time limit by time-sliced evolution sequences. 
Based on that insight, the construction of tunable, noise decoupled, conditional gates composed of 
alternating pulse sequences is proposed. In particular such an approach can lead to a significant 
improvement in fidelity as compared to a strictly periodic control sequence. Moreover we analyze the 
decoherence impact, the relation to the filter for classical noise known from dynamical decoupling 
sequences, and we outline how an alternating sequence can improve spin sensing protocols. 


The interaction of a quantum systems with its envi¬ 
ronment is due both to coherent control fields and en¬ 
vironmental noise. Suppressing the environmental noise 
is crucial for extending the coherence time and for col¬ 
ored noise is frequently performed by dynamical decou¬ 
pling techniques [1-3] . Such approaches aim at refocusing 
undesired couplings by means of a control field or con¬ 
trol pulse sequence, without the need of additional ancil- 
las or measurements. Whereas these concepts are well- 
understood for protecting single isolated qubits, their 
extension to multiple and interacting quantum systems 
turns out to be much more challenging [4-6] . In partic¬ 
ular the question arises how to decouple a qubit from a 
noise background while at the same time retaining very 
specific desired interactions. Such a situation occurs 
in the construction of decoupled quantum gates: Even 
though their feasibility has been theoretically proven [7— 
9], an actual implementation requires concepts that are 
tailored to the specifics of the system. Experimental re¬ 
alizations range from decoupled single [4, 10] and mul¬ 
tiqubit [5, 11, 12] spin gates to single [13] and multi¬ 
qubit gates [14] in trapped ions. Another example is 
given by dipolar recoupling techniques in NMR exper¬ 
iments [15, 16], which aim at the controlled and selective 
reintroduction of specific interactions for structural in¬ 
formation gain. 

More generally, control sequences can be used to con¬ 
struct and design filters [6, 17-20] transmissive exclu¬ 
sively for very specific interactions and frequencies. This 
insight has lead to the construction of quantum ‘lock¬ 
in’ setups [21, 22] or the design of Hamiltonians [23] by 
choosing appropriate control parameters. Moreover de¬ 
coupling sequences can be intuitively interpreted in the 
filter framework [6, 17, 24], with decoherence arising as 
the overlap of a filter with the corresponding noise spec¬ 
trum. Importantly, scanning a narrowband filter in fre¬ 
quency allows for the tomography of the noise spec¬ 
trum [25, 26], the detection of ac-helds [27, 28] and even 
the sensing of (single) individual spins [29-33] . Further¬ 


more, the insight following from filter descriptions has 
been used for the optimization of quantum memories [34] 
and for the construction of sensitive mass [35] and me¬ 
chanical motion [36] spectrometers. 

Here we analyze the filter description for the coher¬ 
ent coupling of hybrid spin systems. As an explicit ex¬ 
ample, though not limited to, we consider the nitrogen 
vacancy center (NV _ -center) in diamond, characterized 
by an electronic spin-1 ground state, and its (hyperfine) 
coupling to individual (nuclear) spins. Both the detec¬ 
tion [29-33] and manipulation [12, 30] of single nuclear 
spins, e.g. 13 C, 29 Si or 14 N nuclear spins, via the cen¬ 
ter spin have been demonstrated. Two relevant regimes 
emerge and are well-described in such a filter formalism: 
the weak and strong coupling limit which correspond to 
close [5] and distant nuclear spins [12, 29-31] from the 
center, respectively. 

Section I is devoted to the filter description in the weak 
and the strong coupling regime, which represents a good 
description for short interaction times and hence small 
rotation angles. This description reveals universal prop¬ 
erties around the resonance points of the filter, which 
allow for tuning of the interaction type and strength. 
Section II then shows how to extend the applicability of 
these concepts beyond the short time limit and reveals 
the ability for tunable gate interactions based on alter¬ 
nating control pulse sequences. An application of this 
insight is provided in section III, linking first the gate for¬ 
malism to the decoherence description, followed by pro¬ 
viding concepts for the use of the alternating sequences 
in spin sensing protocols with improved frequency reso¬ 
lution. Appendix A analyzes in detail the limitations of 
the filter formalism. Appendix B provides the derivation 
of the filter formulas and Appendix C an analysis of the 
strong coupling limit for a spin 1/2 control spin. Last, 
Appendix D reviews the filter-grating description which 
simplifies the analysis and filter calculation for alternat¬ 
ing pulse sequences. 
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I. FILTER FORMALISM FOR COHERENT 

SPIN INTERACTIONS 

We will consider the generic Hamiltonian 

H = 1 a z + j a x ® S z (t) (1) 

describing the interaction of a nuclear (‘target’) spin with 
energy separation w, hyperfine coupled with frequency A 
to an electronic (‘control’) spin. Here Uk (k£{x, y, z}) de¬ 
notes the Pauli matrices acting on the target qubit, and 
for clarity the analogs on the control spin are denoted by 
s k in what follows. S z (t) describes a z-type coupling on 
the control spin, projected on the (pseudo)spin-l /2 sub¬ 
manifold as selected by a sequence of population invert¬ 
ing 7r-pulses; the latter reflected in the time dependence. 
Such a form holds for any type of spin couplings with 
significantly different energy scales which lead to the ex¬ 
clusion of population exchange with the control spin sys¬ 
tem. In particular we assume the control spin to be of 
spin-1/2 or spin-1 as in the nitrogen vacancy (NV) center 
in diamond and the target (nuclear) spin to be a spin- 1/2 
which applies to a wide variety of potential nuclear spin 
targets. It is important to note that the target spin is 
not restricted to spin- 1/2 and higher spins can be intro¬ 
duced straightforwardly into the formalism by replacing 
the spin matrices by their higher order analogs. 

For the control spin two specific configurations are 
considered, dependent on the two-level (sub)-system in¬ 
volved in the control sequence: A spin-1/2 control qubit 
or the spin -1 sub-manifold {m s =±l}, both denoted by 
‘ spin-1/2 control spin ’ in the following and describable by 
S z (t) = s z s(t). A spin -1 sub-system composed of m s = 0 
and one of the states m s = ± 1 , referred to as 1 spin-1 con¬ 
trol' ’, leading to S z (t) = (1/2)(s z s(f) ± 1) with s z and 
1 the Pauli z-matrix and identity defined on the cou¬ 
pled sub-manifold, respectively. Herein the stepfunction 
s(to) = 1 at the initial time with a sign change at each tt- 
pulse time t k , i.e. s(t) = X/^ = 0 (-l) fc ^( 4 +i - t) d(t - t k ) 
for Wpulses with the final time tf = t_/v+i and the Heav¬ 
iside stepfunction 9. 

The effective energy separation of an environmental 
spin with Larrnor frequency coo = 7 B (H$ = (l/2)u;o <x) 
coupled to a control spin -1 by the hyperfine interaction 
A (Hhyp = Act (g) l/2(s z s(t) ± 1) is characterized by 
the dressed states of the unconditional (‘s,-independent’) 
contributions. Thus, uj = |u ; 0 ± A\ [30, 37], which allows 
for distinguishing even target spins of the same species 
by their position dependent coupling properties. The 
orthogonal hyperfine component in ( 1 ) is then just the 
orthogonal component with respect to the dressed en¬ 
ergy levels; parallel components can be neglected in the 
perturbative limit evolution considered below. Adding a 
continuous resonant Rabi driving fl on the target spin as 
has been performed in [5], can be described by the Hamil¬ 
tonian H = (fi/2) (jx + (A/2) o z ® S z (t), and thus corre¬ 
sponds to a global coordinate rotation of the generic form 
(1) by replacing toa z —> flcr X: Aa x —> Aa z and a y —> —a y . 


In the following the hyperfine coupling induced inter¬ 
action is analyzed both in the strong and in the weak 
coupling limit under the influence of a periodic CPMG 
IV-pulse control sequence [5, 38] of pulse separation 2r 
and total time t = 2 Nt, i.e. t k = — r + k (2r) (fce[l,N]). 
The effective coupling strength can then be described by 
a control pulse dependent filter J- with resonances char¬ 
acterized by widths 1/t in frequency. Such a description 
is based on a first order Magnus expansion and holds true 
in the limit of small total evolution angles as shown in 
Appendix A. The extension to much longer timescales is 
discussed in the subsequent section. Importantly, besides 
allowing for specific frequency selective gate interactions, 
the CPMG control sequence leads to a decoherence de¬ 
coupling [3, 39], e.g. extends the coherence time T 2 for a 
quasi-static noise bath to [ 20 ] = N 2 ^T 2 . 


A. Weak coupling limit A <C 10 


The weak coupling limit is essentially insensitive to the 
control spin being either of the ‘spin- 1 / 2 ’ or ‘spin- 1 ’ type, 
the latter requiring simply a substitution A —» A/2 in the 
subsequent expressions as a result of its effectively weaker 
hyperfine coupling. The total conditional evolution in a 
rotating frame with respect to H 0 = ( co/2)a z follows as 

Hint = exp 1 N ) ® s z + 0 ([Af] 2 )J ( 2 ) 


with the frequency selective filter given by (see Ap¬ 
pendix B) 


d~w [ojt , N) 


2 

N(ojt) 


sin 2 (wr/ 2 ) 

cos(o;t) 


sin(A r wr + { 7 r/ 2 }) 


( 3 ) 

and the argument in curly brackets being absent/present 
for N even and odd, respectively. Such a filter defines the 
effective coupling frequency as (A T w ) and is illustrated 
in figure 1 (a). It reveals sharp resonances around lot = 
(2 k + 1) 7 r /2 (k £ INg) with a frequency width inversely 
to the total time, and a total peak height decreasing as 
oc 1 /{lot). The rotation axis a$ in the resonance region 


lot = (2k + 1) 7t/2 + ctt/(2N) (4) 


follows a simple linear behavior, that is 
<j 0 = (—l) fc [cos 4>a x — sin (f><J y \ with <j> — c-{n/2 ), (5) 


where c € (—2,2), except for iV=l & fc=0 in which case 
the restriction c £ (—1,1) holds. This angle dependence 
is valid in any limit of N. Moreover in the limit of large 
N (N >4), or for any N at c=0, the filter (3) takes on 
a universal behavior around the resonances 


•d~ w (d ^0 


4 | sin(c 7 t/ 2 ) | 
7 r 2 |c| (2k + 1 ) 


( 6 ) 


resulting in the filter amplitudes F w (0, k) = (2/ 7 r) 1/(2 k+ 
1 ) for a c =0 conditional a x , and F w (l,k) — 
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(4/7T 2 ) 1/(2 k + 1) for a c=l cr^-type rotation. Impor¬ 
tantly, the N independence in (6) ensures a linear ac¬ 
cumulation of the rotation in time. On the other hand, 
the k dependence reflects a reduced coupling for higher 
order resonances. In the validity range of (6), the total 
rotation angle for a given c and N turns out to be in¬ 
dependent of k, however associated with a longer total 
time for increasing k. Moreover, the filter is symmetric 
in c exclusively in the limit (6) of large N. 


B. Strong coupling limit A uj 

In the strong coupling limit the achievable interactions 
crucially depend on the spin properties of the control 
qubit. Only for the ‘spin-1 configuration’ a conditional 
interaction that is scaling linearly in time can be achieved 
and will be focused on in the following; the ‘spin-1/2 
configuration’ leads at most to a scalable (less relevant) 
unconditional evolution and is briefly outlined in Ap¬ 
pendix C. 

For a ‘spin-1’ control qubit the evolution in a rotating 
frame with respect to the hyperfine Hamiltonian can be 
described by 


U int (t) = exp (-i [F c (A t, N) a £ (8) s z 

+F u (Ar,N)a%® 1] +0(M 2 )) 


(7) 


connected to the the original (1) frame evolution by 
U(t n ) = exp(=F?'(A/4)f„ cr x ,) iTint(in) at the end of a ba¬ 
sic decoupling sequence t n = to + 2 nr. Upper and lower 
signs here and in the following are in accordance with 
S z (t) = (1/2) (s z s(t) ± 1), dependent on the spin-1 sub¬ 
system connected by the control pulses as previously de¬ 
fined. Thus the evolution is characterized by both con¬ 
ditional and unconditional contributions, with the cor¬ 
responding filters of width A A a 1/fa 1/N defined as 
(see Appendix B) 


•7c (At, TV) 
•7/(At, TV) 





FIG. 1. Gate filter for periodic control pulses. Filter func¬ 
tion T describing the coherent spin interaction in the limit of 
small total evolution angles and for N = 6 periodic 7T-CPMG- 
pulses with interpulse delay 2 t on the control spin, (a) Weak 
coupling limit A <C u, which leads to a purely conditional 
interaction as follows out of (2). The corresponding rotation 
axis cr^> around the resonance peaks as given by (5) is illus¬ 
trated in the inset. Amplitudes for higher order resonances 
decay as 1/(wt) (black dash-dotted), whereas the resonance 
peak width is inversely proportional to the total time, (b) 
Strong coupling limit d > u with both conditional (blue) 
and unconditional (green) resonances and filters. For both 
(a) and (b) red dashed lines denote the universal amplitude 
behavior (6) and (11) valid in the limit of large TV’s. 


with <j> = c7t/2 for c £ (—2,2), except c £ (—1,1) for 
and N= 1, and restricted k dependent c-ranges for cr^ 
and N = 1,2. 

For large N (N > 4) both the conditional and uncondi¬ 
tional filter are well-approximated by the universal form 


and illustrated in figure 1 (b). Again the contribution in 
wavy brackets is present exclusively for N being odd. 

The resonance (region) condition, in a notation analo¬ 
gous to (4) is given by 

(At) c = (2k + l)n + CTr/N, 

(At) u = k (2n) + cn/N 

for the conditional (c) and unconditional (u) filter, re¬ 
spectively. The corresponding rotation axis follows as 


cr% = =F (cos(</>) ct 2 ± sin(</>) a y ) 
cr)) = cos (<j>) cr~ ± sin <jxr y 


( 10 ) 


•7) 


c/u 


sin (c7r/2) 

C 7T 


(ii) 


which exactly holds for any N on the resonance c=0. 
Thus, c=0 leads to a ov-type rotation with T c / U = 1/2, 
whereas |c|=l corresponds to a u y rotation with T c / U — 
1/7T. The N independence of (11) indicates the additivity 
in time; the additional k independence reveals that the 
coupling amplitude, despite a modification of the ability 
for noise decoupling, is preserved for higher order res¬ 
onances. Out of (8) and (11) it directly follows that 
the conditional filter is symmetric in c in any limit of 
A, whereas this holds true for the unconditional filter 
exclusively for sufficiently large N. The (un) conditional 
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FIG. 2. Sliced evolution and long time limit extension (weak coupling limit A <C cj). (a) Effective coordinate system rotation 
by ft — C7r after a single slice with even pulse number n, tuned to the resonance parameter c (4) and characterized by an angle 
tj> = C7T /2 of the rotation axis (5). (b) Interpretation of a large N sequence on resonance c beyond the filter validity by 
decomposition into s n'-pulse slices. Each slice is followed by a coordinate rotation ft' = c'n with c' = an!/N up to a total 
angle ft = C7r. Thus the rotation axis (red arrow), characterized by <j> = c'(n/2) in the local frame, is globally rotated in time, 
(c) Numerical simulation of a IV=20 7r-pulse sequence for A/uj = 0.06, resonance order fc=0, and the initial state | J. +) with 
the control spin states |±) = 1/V2 (| 4) ± | t))- Solid lines follow out of 10 alternating ±c (c = 1) slices as illustrated in (d) 
with rij= 2 7r-pulses each. Dashed lines describe an equidistant (+c) CPMG sequence of 20 7r-pulses. The total time shown 
corresponds to an expected 7r/2-pulse evolution based on the filter formalism. The difference in total time follows from the 
filter asymmetry in c in the limit of small pulse numbers n. Red arrows indicate the 7r-pulse positions for the sliced (solid line) 
evolution; the different slices are shaded in different colors in analogy with the sequence coloring in (d). The final fidelity for the 
conditional n/2 evolution can be extracted as F = 99.9% (±c sliced) and F = 90% (CPMG). (d) Alternating ±c control-pulse 
sequence composed of s slices. 


contribution at the unconditional (conditional) resonance 
peak does not scale in N and is thus negligible for a suf¬ 
ficiently large number of pulses; moreover it is zero for 
c=0 and |c|=l for an even and odd number of pulses N, 
respectively, in any limit except for AT=1[40]. 


II. SLICED GATE EVOLUTIONS AND 
ALTERNATING PULSE SEQUENCES 

The filter description and its associated properties dis¬ 
cussed previously hold true as long as the interaction 
time and thus the evolution angle is small 0 <C 1 (see 
Appendix A). For gate interactions this forms a severe 
limitation, which can be overcome by alternating pulse 
sequences as shown in the following. 

We will assume that the total evolution can be split 
into s slices [tj,tj + 1 ) (j £ [0, s-1]), with slice j consisting 
of rij-periodic CPMG-control pulses with characteristic 
timescale r ? as shown in figure 2 (d). Moreover i 0 =0, the 
total time t=t s and the slice time Atj=tj+i~tj. Then the 


total evolution can be cast into the form 

s —1 s-1 

U = e~ l H Atj = e -1 H ° * ]^[ e* H ° tj U^2 (Af^ )e _I H ° tj 

3=0 3=0 

( 12 ) 

where Hq = (w/2) a z and Hq = ±(A/4) a x for the weak 
and strong coupling configuration, respectively. U^t de¬ 
notes the rotating frame evolution operator as defined 
in (2) and (7) with the index (j) indicating that the re¬ 
placement s z — y — s z has to be made in these expressions 
for an odd number of preceding pulses. Assuming the 
rotation angle per slice to be small, its evolution is well- 
described in the filter description; e.g. for A<Cw the rota¬ 
tion angle 6j=2(A/uj) (rijUiTj) !F w (ojTj , rij), which is small 
as long as rij is. Following (12) the only modification 
consists of an effective coordinate exp(iHaAtj) rotation 
between subsequent slices j and j + 1. 

The rotation of the coordinate frame is closely linked 
to the rotation axis angle of the evolution. This property 
is analyzed for the weak coupling limit in the following, 
however the same conclusions hold true for the strong 
coupling counterpart. After a single slice j around a ro¬ 
tation axis = cn/2 following out of conditions (4) and 
(5), the local coordinate system can be interpreted as 
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being rotated by 

1 9j = 2rij (i 'jJTj ) = (2k + 1 )rijTr + err (13) 

for the subsequent slice. This rotation, for rij being even, 
is just twice the rotation axis angle §j = 2<j> (see fig¬ 
ure 2 (a)). An additional 7r-inversion occurs for rij be¬ 
ing odd, i.e. r)j = 2(j> + n. However in that latter case 
also s z — } — s z as outlined previously by the index (j) to 
the evolution operator, and thus the rotation axis is ir- 
inverted simultaneously and the subsequent conclusions 
remain unchanged [41]. 

Thus, to globally maintain a fixed rotation axis </>, the 
subsequent slice has to be tuned to — c. An alternating 
evolution sequence c —► —c — >■ c —► — c... as depicted in 
figure 2 (d) therefore allows for significant rotation angles 
beyond the perturbative limit around an arbitrary but 
fixed rotation axis (5) as characterized by <p = c7r/2. This 
is illustrated for a 7r/2 evolution in figure 2(c), compar¬ 
ing a non-alternating N=2Q pulse sequence for c=l to 
a c = ±l alternating sequence of s=10 slices with rij =2 
pulses each. Whereas the former simulation shows a clear 
deviation from the filter based expectation as a result 
of the small-angle approximation break-down, the latter 
sliced evolution leads to an (almost) perfect cx y -type con¬ 
ditional 7r/2-rotation with fidelity F = 99.9%. 

For an even number of alternating ±c slices (or for 
c=0), exp(iHot) oc 1 such that the rotating frame evo¬ 
lution is equal to the original one. Moreover in this case 
and for rij=n, the effective ‘filter’ for the c-axis rotation 
in analogy to (2) is given by 

Jos = 1/(2n) ((c + n) J +c + ( c + n) J_ c ), (14) 

or simply J e ff = J~± c whenever J is symmetric in c. 

Remarkably the c=0 resonance peak is purely additive 
even in the long time limit without any adjustments, a 
fact that can already be seen in the ^-independence of 
the corresponding resonance condition (4). On the con¬ 
trary, two subsequent c=+l slices would be subtractive 
to first order as the coordinate frame rotation inverts the 
global rotation axis for the second slice. 

The above insight can also be applied for the interpre¬ 
tation of the long time (large N) evolution of equidis¬ 
tant CPMG sequences (t, = r). For that purpose, an 
.ZV-pulse sequence characterized by a resonance parame¬ 
ter c according to (4) can be formally sliced into s seg¬ 
ments of n' = N/s pulses each. The resonance and ro¬ 
tation properties for each slice are then characterized by 
d = (n'/N)c, as follows by comparing the resonance con¬ 
dition per slice out = (2k + 1 )t/ 2 + d-rr/(2n') to (4). 
Choosing sufficiently many slices such that the rotation 
per slice is small allows for the interpretation within the 
filter formalism. Namely, as illustrated in figure 2 (b), the 
rotation of each slice with d is followed by a coordinate 
system rotation J = d 7r; up to a total rotation {) = ck 
for the total sequence. Essentially this equals a rotation 
of the rotation axis in time. Therefore, based on the slic¬ 
ing approach, the filter description retains its usefulness 
beyond the perturbative limit. 


III. SPIN SENSING AND FREQUENCY FILTER 
FOR ALTERNATING PULSE SEQUENCES 

The conditional coherent interaction induced by 
Hamiltonian (1) leads to a mutual, time-periodic entan¬ 
glement of the control and target spin. This effect can 
be observed as a coherence decay on the control spin 
subsystem, which along with the frequency selectivity of 
the control pulse sequence is widely used for sensing in¬ 
dividual environmental spins [29-31]. More precise, the 
coherence follows as (with s + the ladder operator for the 
control spin) (s+(t)) = (s+(0)) tr[^ =+1)t ^ ( n s t 5= ^ 1) p c ], 
which in the A <C cu limit and assuming the target initial 
state p c = (1/2)1 to be completely mixed, leads to 

(s + (t)) ~ (s+(0)) cos (AtF w (ouT, N)) . (15) 

For sufficiently short times, the coherence evolution fol¬ 
lows as (s+(f)) ~ (s+(0)) exp(— (1/2) A 2 t 2 F w (uiT, N) 2 ). 
This latter expression is formally equivalent to the decay 
under classical noise [29] , with the corresponding filter 
Jclass = Jj [17, 20] and assuming a noise bath with dis¬ 
crete spectrum S(ou') = 6(ou — ou'). 

The validity range of (15) is limited by the applicability 
range of the filter formalism. In practice, as simulated for 
a CPMG sequence in figure 3 (a), the coherence is well- 
described by (15) within the initial sign inverting regime 
(s+(0 ))—>— (s+(0)), i.e. up to a rotation angle n of the 
conditional evolution. It forms the typical regime used 
for spin-sensing. Going beyond this timescale reduces 
the dip depth of the main resonance again, whereas the 
vicinity regions still increase in depth. This then leads to 
an oscillatory broadened regime associated with a loss of 
frequency resolution. More precise, such a regime makes 
it hard to distinguish different spins of similar frequency 
as a result of their overlapping signal envelopes; in addi¬ 
tion to the drawbacks associated with the (potentially) 
reduced coherence dip amplitude and the more intricate 
signal pattern. Moreover the filter formula description 
(15) breaks down at these timescales; except for c=0 
(wto=7t/2 in figure 3 (a)) as is clear out of the ‘slicing’ 
analysis of the preceding section. Thus there exists a 
trade-off between an anticipated increasing sensing res¬ 
olution in time out of the filter description (Aw oc 1/t) 
and an upper limitation by the onset of the oscillatory 
broadening at longer times; besides the limitation on the 
coherence time set by the classical noise background. 

Extending the time until the onset of the oscillatory 
regime can be achieved by lowering the effective coupling 
amplitude. That way the filter description and its time 
inverse frequency width scaling retains its validity on a 
longer timescale, which, provided the (decoupled) coher¬ 
ence time exceeds that timescale, enables to further in¬ 
crease the sensing resolution [37]. Thus, for a fixed total 
time, appropriately tuning the coupling strength allows 
to turn a signal with previously oscillatory broadened 
envelope and reduced dip depth, into a clearly distinct 
peak structure of larger or maximal peak depth, each 
peak width scaling with the inverse total time. 
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FIG. 3. Spin sensing (weak coupling A/uj = 0.05). (a) Control spin coherence (s+) (red) for a CPMG pulse number IV = 24 
(upper) and IV = 48 (lower) with total time t = 2 Ntq. Black lines indicate the coherence evolution as expected out of the filter 
formalism (15). The squared filter T 2 is illustrated for comparison, (b) Coherence for an alternating c = ±1 IV = 48 pulse 
sequence with M= 12, 2n=4 (red) and M = 6, 2n=8 (blue, dashed). Black lines correspond to the filter expectation for the 
M=12 case. Lower panel: Filter (17) for the M=12 (left) and M=6 (right) configurations and its composition (in arbitrary 
units) out of grating (dotted) and block filter (dash-dotted). 


Such an amplitude damping effect is realized by higher 
order resonances (fc>0) [30], by superimposing a contin¬ 
uous Rabi driving on the control spin [42] , or by lowering 
the effective coupling amplitude by gradually reducing 
the periodicity of the control sequence [37] . This latter 
concept can be achieved by an IV-pulse sequence com¬ 
posed of s alternating ±c slices with n pulses each (fig¬ 
ure 2(d)), choosing the inter-pulse timescales such that 


T±c =T 0 ±C 


TO 

(2k + 1) n 


(16) 


Here r 0 = (1/2) (r +c + r_ c ) or simply r 0 = t/(2N) for 
an even number of slices s. That way, at the resonance 
condition ujtq = (2k + 1) n/2, the system is alternatingly 
driven at the ±c resonance, i.e. uit± c = (2k + 1)tt/2 ± 
cn/(2n). This corresponds to an (approximate) single 
axis rotation with a well-controlled and reduced coupling 
amplitude as discussed in the preceding section. 

However, a reduced periodicity comes at a price, 
namely additional peaks appearing in the sensing spec¬ 
tra. This is best understood by noting that the filter for 
an M times periodic repetition of an m-pulse block can 
be described as the overlap of a grating Q with a block 
filter J^iock [23, 37] (see Appendix D) 

T±c(uf = ^ G(ujt, M).Fblock(wTo,?Tl,c) 2 . (17) 

We will assume the number of slices s to be even, in 
which case M = s/2 and m = 2n with a ‘block’ consist¬ 
ing of a single ‘c —>• — c’ sequence. The grating is then 
given by (m even) Q = sin 2 (wI/2)/sin 2 (o;t/(2M)) and 
the block filter is calculated analogue to the CPMG filter 


(see Appendix B) by just adapting s(t) to the modified 
pulse sequence. Thus the total filter can be interpreted as 
the overlap of a grating with peak width inversely to the 
total time A uj oc 1 /t modulated by a broader block filter 
of width Aw oc (t/M)~ x . Compared to a perfectly ‘sym¬ 
metric’ grating ( M=N , ro=1), in which the grating con¬ 
sists of peaks wr 0 = (2fc + l)7r/2, (m-1) additional peaks 
per resonance order k and separated by A(wr 0 ) = jr/m 
emerge for M = N/m, gradually introduced into the fil¬ 
ter by reducing the symmetry of the ?n-pulse blocks [37] . 

Figure 3 (b) illustrates the filters for a sequence of ±c 
alternating slices, a fixed total pulse number IV=24 and 
different pulse numbers n per slice. Previously located in 
the oscillatory regime for a fully periodic IV=24 sequence 
(figure 3 (a)), the alternation leads to well-defined peaks 
of frequency width given by the inverse total time and 
well-described in the filter framework; though at the ex¬ 
pense of the emergence of additional peaks for a given 
frequency. 


IV. CONCLUSION 

We have shown how the evolution operator of interact¬ 
ing spins subject to a control pulse sequence can be de¬ 
scribed in a filter description. This holds true for (hyper- 
fine) coupled spins without population exchange in both 
the limiting cases of strong and weak coupling. A uni¬ 
versal behavior can be ascribed to the resonance region, 
allowing for a straightforward identification of both the 
interaction type and amplitude. Importantly, it has been 
shown that the interpretation in the filter framework can 
be extended to the previously inaccessible limit of long 
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times and large rotation angles. Thus its insight can be 
used for the construction of tunable decoupled quantum 
gates. It should be noted that obtaining the exact solu¬ 
tion of the conditional spin evolution is a straightforward 
task [5, 29-31], though its formulation is more involved, 
and quickly becomes a numerical task involving the in¬ 
version of trigonometric functions. In contrast, the fil¬ 
ter formulation provides an intuitive approach and com¬ 
bining both methods can lead to further optimizations 
whenever the limiting cases of strong and weak coupling 
are not strictly fulfilled. Moreover, we have analyzed the 
impact on the (control spin) subsystem coherence. Re¬ 
markably, the sensing relevant regime is well-described in 
the filter formulation and is closely related to the filter for 
classical noise. Furthermore it has been shown that sens¬ 
ing resolution can be improved by alternating sequences 
with well defined coupling amplitudes as follow directly 
out of the filter description. 


ACKNOWLEDGMENTS 

This work was supported by an Alexander von Hum¬ 
boldt Professorship, the ERC Synergy grant BioQ and 
the EU projects DIADEMS, SIQS and EQUAM. 


Appendix A: Limitations to the filter description 


Hamiltonian (1) in a rotating frame with respect to 
H 0 = (oj/2)a z , well suited for the filter analysis in 
the weak coupling regime A <C w, follows as H lnt = 
(A/2) exp (iojt) + h.c.] 0 S z (t). The evolution oper¬ 
ator in the Magnus expansion can then be written as 


C/int = exp 



(Al) 


with the first two expansion contributions 


Hi(f) 


—i A/2 [ H int (t')dt' 
Jo 



Hint (£2) 

(A2) 


and the definition H lnt = (A/2) H lnt . A completely ana¬ 
logue treatment can be carried out for the strong coupling 
limit in a rotating frame with respect to the hyperhne 
contribution. 

A sufficient condition ensuring absolute convergence 
of the sum appearing in (Al) is given by [43] 
/o ||-H r int (t , )||d£ / < 7 r, with the matrix norm here and in 
the following defined as the operator norm. 

As the filter description is based on the first order con¬ 
tribution alone, its validity is restricted to timescales on 
which higher order contributions are sufficiently small. 


With ||Ffi n t (t) || = J, a time independent quantity as a 
result of the unitarily invariant matrix norm, it follows 
that ||fli(t)|| < (A/2)Jt and ||n 2 (f)|| < 1/2 (A/2) 2 t 2 J 2 
and under the assumption of absolute convergence [7] 


E n *(*) 

i=k 



(A3) 


with a constant Ck = 0(1). 

Recalling the filter definition (see Appendix B) F = 
(l/t)\ f 0 exp(iujt')s(t')dt'\, the first order contribution 
takes the form ||f2i(£)|| = \ A\/2JFt and the above esti¬ 
mations allow for an upper bound on time for the filter 
description validity. 

More approximate, for sufficiently many pulses N and 
close to resonance 


mm < (4 


2 ,t 


dfi 


H ln t(t\ 


dfg-Hint (t 2 ) 


<('!) J 2 F [dt^ = 


i o 

A \ 2 1 
— Jt 1 -F 
2 2 


(A4) 


where the ‘coarse grain’ estimation WJ^ 1 dt 2 H int (t 2 )\\ < 
FJt i, approximately valid in the quasi-resonant regime 
for sufficiently long times (large N), has been used. Thus, 
||H 2 (f)|| ((A/2) Ft J) 2 , which compared to the first 

order contribution and its associated rotation angle 6 
||fli(t)| = (A/2FtJ) = 0/2 J, reveals that the second 
order contribution is of O([0/ 2] 2 ). Thus, the general va¬ 
lidity of the filter description is limited to small rotation 
angles 0; except for the special cases discussed in the 
main text (e.g. c= 0 or alternating sequences). 

In the weak coupling limit, the evolution operator (2) 
up to second order takes the form I7i nt ~ exp(— i[0/2 
s z + (0/ 2) 2 (F^/F 2 ) a z 0 1]) with 0 = AF w t the first 
order rotation angle, and F w and F^ the first and sec¬ 
ond order filter, respectively. These filter contributions 
are compared in figure 4, verifying that in the resonance 
region F^/F w < 1 and thus indeed the second order 
contribution < (0/2) 2 as expected out of the previous 
discussion. 



FIG. 4. Comparison of the first and second order filter 
contribution in the weak-coupling limit for a number N= 6 of 
periodic pulses. 























Appendix B: Filter function derivation 
1. Weak coupling limit 


The rotating frame evolution, retaining only the first 
order Magnus expansion contribution, takes the form 

C/ in t ~ exp^-i^ [cr + J e lult s(f , )df , +h.c.](g>s, z ^ (Bl) 

with s(t)G{ — 1,1} the control pulse jump function as de¬ 
fined in the main text. Thus it essentially remains to 
calculate x(t) = f* exp(iujt')s(t')dt' and the filter (3) 
and rotation angle (5) subsequently follow as T w = 
(1 /t) |x(t)| and tan</> = Im(x(t))/Re(x(t)), respectively. 

Using that the pulse times for an equidistant TV-pulse 
sequence are given by t & = — r + k (2r), with start and 
end time to = 0 and fjv+i =t = 2 Nt 

N 

x(i)=E(- 1 ) fe 

k—0 

(- 1 ) + (_i)^ e “ (2ArT) - 2e- iwT E [-e iui2T ] k ) 

k =1 / 

(B2) 




and upon evaluation of the sum using the geometric series 
property 


x(t) 


i 

uj 




gicuATr 

cos(wr) 


f —isin(u;A f T)'l \ 
\ cos (ojNt) J ) 


(B3) 

with the upper and lower case corresponding to N being 
even and odd, respectively. 


2. Strong coupling limit 

Hamiltonian (1) in the rotating frame with respect to 
the hyperfine contribution and for the ‘spin-1’ control 
spin configuration S z (t) = 1/2 (s z (t) ±1) (corresponding 
to upper and lower signs in what follows) takes the form 

H int = | (s+e 4 ^ [s * sRt)±t] + h.c.) (B4) 

with c+ the ladder operator in the cr x eigenbasis, i.e. <;+ + 
h.c. = a z and i + h.c. = a y , and s F (f) = f 0 s(t')dt'. 
For a periodic (CPMG) iV-pulse sequence s F (f) can be 
expressed as 


A separation of (B4) into a conditional and unconditional 
part is obtained by noting that 


cos([A/2]s F (f)s z ) = cos([A/2]s F (t))l 
sm([A/2]s F (t)s z ) = sin([A/2]s F (f)) s z . 


(B6) 


The total evolution in the rotating frame and approxi¬ 
mated by the first order Magnus expansion then follows 
as 


Hint — exp ( - i — 


with 


(?+ X* + h.c.) ® s z + (?+Xi/ + h.c.) (g) 1 


(B7) 


Xc = * J o dt'sin^ flF (0) e ±t(A/2)t ' 

Xt= jfVcos^MO)e ±i(A/2)t '. 


(B8) 


The corresponding conditional and unconditional filter as 
defined in (7) and (8) then follow as J- c / u = (l/t)|x//J = 
(l/t)|x~/J) and the rotation angle as defined in (10) is 
given by tan</> = Im(x+/J/ Re (x(!/J’ 

Thus it remains to evaluate the expressions (B8). In¬ 
serting the expression (B5) for s F (f) 


*Ms- i >7’‘ +V( 

2 k =0 J *k 


i(A/2)(t'-k(2 T )) 


- h.c.)e 


±i(A/2)t' 


(B9) 

and analogue for x« • Performing the integration, insert¬ 
ing the expressions for t fc =—r + k( 2r), separating the 
‘special times’ to & Uv+i = t, and finally using the ge¬ 
ometric series property in close similarity to the weak 
coupling treatment (B2), results in the expressions 


xt = [1 ^ (—1) N e ±iANT ] 


l T 

2A T 2 


0 ±i(A/2)Nt 


l 

T T 


A cos((A/2)r) 


0 ±i(A/2)r 


±i sm((A/2)Nr) 

— cos((_A/2)iW) 

(BIO) 


and 




A sin((A/2)r) 


0 ±i(A/2)r 


sin ( —Nt 


(BH) 


Upper and lower cases in wavy brackets correspond to the 
total pulse number N being even and odd, respectively. 


N 

s F {t) = E(-!) fc ~ k 2r ) 0 (^+i - t) #(* - *fc) (B5) 

k =0 


Appendix C: Spin-1/2 control spin in the strong 
coupling limit 


with the pulse times tf. = —r + k (2 r) and the initial and For the strong coupling limit i»w and assuming the 
final time fo = 0 and fjv+i=t, respectively. control qubit to be of the ‘spin-1/2 type’, the rotating 
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frame Hamiltonian is given by 


coupling limit is given by 


H int = | (<t+e lAs -« s *+h.c.) (Cl) 


Ht) 


1 

t 



e iut ' s(t')dt' 


(Dl) 


with and sp’(t) = / Q f s(t')dt’ as defined in Ap¬ 
pendix B 2. 

Thus, the rotating frame evolution follows as 

Hint = exp t \T c a y <8 s z + T u a z ® 1] + £>((wt) 2 )^ 

(C2) 

with T c = (l/i)Im(£), T u = (l/f)Re(C) and the defini¬ 
tion ( = f* dt' exp(iAsF(t')). With (B5) and using the 
same procedure as in Appendix B the filters are readily 
evaluated to 


{ 1/ (N At ) (l — cos(At)) } (C3) 

T u = 1/ (At) sin (At ). 

Upper and lower cases in curly brackets correspond to 
an even and odd number of total pulses N, respectively. 
Remarkably, the conditional filter T c does not lead to a 
scalable evolution, i.e. it scales inversely with the pulse 
number or is even zero for an even number of pulses. 
More precise, at most a single pulse unit contributes to 
the evolution and thus such a configuration is not suit¬ 
able for conditional gate interactions. In contrast, the 
unconditional contribution is scalable and does reveal a 
discrete peak structure. 


Appendix D: Grating block-filter construction in the 
weak coupling limit 

The filter in accordance with the definition (2) as fol¬ 
lows from a first order Magnus expansion in the weak 


with the step function s(t) £ {±1} as previously de¬ 
fined describing a pulse sequence consisting of M peri¬ 
odic repetitions of m-pulse blocks. Taking into account 
that the time integration involved can be split into M 
periodic time intervals [tk,tk+ 1] with ti = 0, £m+i = t 
and At = —tk= t/M , allows to rewrite (Dl) as 


Tit) 



e iut ’ s{t') d t' 


M 

nt/M 

e iw(k-l) (t/M) 

/ e iult ' sit 1 ) d t 1 

k =1 

10 


(D2) 


where in the last step the periodicity property s(tk+t r ) = 
[(—l) m ] fc+1 s(f') has been used. Out of (D2) and (17) it 
is straightforward to identify the grating as 


t/(wf, M, to) 


Et(-in fc+1 


k =1 


2 


e iu(k- 1 ) (t/M) 


(D3) 


that upon evaluation of the sum using the geometric se¬ 
ries property results in 


G 


sin 2 (aj i/2) 
sin 2 (oji/(2M)) 


cos 2 (ujt /(2 M)) 


sin 2 (u;f/2) 

cos 2 (wf/2) 


m even 

to odd, M even 
to odd, M odd. 

(D4) 


Correspondingly the ‘block filter’ is identified as .Fbiock = 
| J^ M exp^wf') s(t')dt'\. For a ±c block with n pulses per 
slice, i.e. to = 2 n, the block filter follows as ^bio C k(^ n ) = 
| Xn{r+ C ,t = 2 nr +c ) + (-1)” exp(i2nwT +c )x„(T_ c ,i = 
2nT- c )\ with x as defined in(B3). 
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